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1-^ ' Abstract 



A nonrelativistic boundstate formalism used in contemporary calculations is investi- 
gated. It is known that the effective Hamiltonian of the boundstate system depends on 
5^ \ the choice of gauge. We obtain the transformation charge Q of the Hamiltonian for an ar- 

bitrary infinitesimal change of gauge, by which gauge independence of the mass spectrum 
and gauge dependences of the boundstate wave functions are dictated. We give formal 
arguments based on the BRST symmetry supplemented by power countings of Coulomb 
singularities of diagrams. For illustration: (l)we calculate Q up to ©(l/c), (2)we ex- 
amine gauge dependences of diagrams for a decay of a qq boundstate up to 0{l/c) and 
show that cumbersome gauge cancellations can be circumvented by directly calculating 
Q. As an application we point out that the present calculations of top quark momentum 
distribution in the it threshold region are gauge dependent. We also show possibilities for 
incorrect calculations of physical quantities of boundstates when the on-shell matching 
procedure is employed. We give a proof of a justification for the use of the equation of 
motion to simplify the form of a local NRQCD Lagrangian. The formalism developed in 
this work will provide useful cross checks in computations involving NRQED/NRQCD 
boundstates. 



1 Introduction 

Recently there has been much progress in our theoretical understanding of nonrelativistic 
QED and QCD (NRQED and NRQCD) boundstates such as positronium, T and remnant of 
toponium boundstates. Following the idea of nonrelativistic effective field theory proposed by 
Caswell and Lepage |jl|, formalisms necessary for precise descriptions of these boundstates have 
been developed significantly p|-[0. At the same time there appeared many new calculations 



of higher order corrections to physical quantities of both the NRQED [|T2l-[|T4| and NRQCD 



boundstates |T5|-||2^ (boundstate mass, decay width, production and decay cross sections, 
etc.). Despite these developments, a completely systematic formulation necessary for compu- 
tations of these physical quantities in perturbative expansions has not been established yet. 
Among current best technologies the asymptotic expansion of Feynman diagrams |^ seems to 
be most suited for calculations of the physical quantities, in particular for fixed-order calcula- 
tions. In addition, effective field theories are powerful tools in order to sum up various large 
logarithms originating from the widely separated scales inherent in the problems. Efficiencies 



and correctness of various effective theories are, however, still under debate. (See Refs. [|, [10 
for discussions on the current status of the formalisms.) 

A notable characteristic in these new developments is that the conventional Bethe-Salpeter 
equation is no longer being used to calculate the spectrum and wave functions of boundstates. 
Instead, one starts from the non-relativistic Schrodinger equation (of quantum mechanics) with 
the Coulomb potential. Then one adds to the nonrelativistic Hamiltonian relativistic correc- 
tions and radiative corrections as perturbations to obtain an effective Hamiltonian (quantum 
mechanical operator) valid up to a necessary order of perturbative expansion. Effective Hamil- 
tonians used in these new formalisms are known to be dependent on the choice of gauge. 

The purpose of this paper is to investigate gauge dependence of a nonrelativistic boundstate 
formalism used in contemporary calculations 0],||14|,|[T^-[^,|Q, which is also closely tied to 



the potential-NRQCD formalism [0]. Our motivations for the study are: (I) In the present 
frontier calculations of higher order corrections to physical quantities of boundstates, often 
the Feynman gauge is used to calculate typically ultraviolet radiative corrections whereas the 
Coulomb gauge is used to calculate corrections originating typically from infrared regions. 
Although much care has been taken to calculate consistently in each gauge separately gauge 
independent subsets of the corrections, it is desirable to clarify gauge dependences of the 
formalisms actually used in these calculations. (II) We would like to find transformations of 
boundstate wave functions when we change the gauge-fixing condition. We may apply these 
transformations to study various amplitudes involving boundstates. Since a physical amplitude 
is gauge independent, once we know how the wave function transforms, we know how other 
parts of the amplitude should transform to cancel gauge dependence as a whole. This would 
provide a useful cross check for identifying all the contributions that have to be taken into 
account at a given order of perturbative expansion. 

Already some time ago, gauge independence of the mass spectrum of the NRQED bound- 
states was shown and studied in detail based on the Bethe-Salpeter formalism pSf-p^: Ref. |25 



gave a brief discussion; Ref. p6| examined a Feynman gauge calculation of the boundstate spec- 



trum at next-to-leading order very closely and showed that an infinite number of two-particle 



irreducible diagrams contribute in this gauge, which in the end all cancel due to fairly com- 
plicated gauge cancellations (This feature is much more complicated in comparison to the 
calculation in Coulomb gauge.); Refs. [^ gave formal arguments as well as perturbative anal- 
yses which apply to all orders of perturbation series. 

In comparison to these earlier works, new achievements of the present work are: 

• We use the BRST symmetry to formulate our arguments, which allows us to treat both 
the NRQED and NRQCD boundstates on an equal footing. In particular we are able 
to discuss gauge dependence of the NRQCD boundstate formalism rigorously using this 
symmetry. 

• Presently, there exist several different definitions of an effective Hamiltonian beyond 
leading order. We introduce an effective Hamiltonian defined naturally in the context 
of time-ordered (or "old-fashioned") perturbation theory of QED/QCD. Then we ob- 
tain a transformation charge Q (quantum mechanical operator) such that the effective 
Hamiltonian and the boundstate wave function change as 

8H,s{P') = [H,s{P') - P°] iQ{P') - iQ\P') [H^siP') - P'] , 
6ip = —iQ-ip 

when the gauge-fixing condition is varied infinitesimally.Q Also, gauge independence 
of the spectrum is shown using the transformation. We define the charge Q of the 
effective Hamiltonian directly in terms of the BRST charge and the field operators in the 
QED/QCD Lagrangian. 

• For illustration: (l)we calculate the transformation charge Q at next-to-leading order; 
(2)we demonstrate gauge cancellations among diagrams by examining an infinitesimal 
gauge transformation of the amplitude for a qq boundstate decaying into q'q"W^W~ . 
From the latter example, one can deduce that the present calculations of the top mo- 
mentum distribution in the ti threshold region at next-to-next-to-leading order are gauge 
dependent. 

Another subject of this paper is to study the problems in a determination of the effective 
Hamiltonian from the on-shell scattering amplitude of a fermion and an antifermion. Generally 
a fermion and an antifermion inside a boundstate are off-shell so that use of a Hamiltonian 
determined in the on-shell matching procedure may lead to incorrect calculations of the physical 
quantities of the boundstate. We clarify this point. 

The same problems do not occur if we use a local NRQED/NRQCD Lagrangian and deter- 
mine (Wilson) coefficients in the Lagrangian by matching onto the full theory, i.e. by matching 
on-shell amplitudes of all relevant physical processes to those of perturbative QED/QCD in 
nonrelativistic regions. In this case one should calculate amplitudes for a number of processes 



* One may conjecture that gauge dependence can be described by unitary transformation, since effective 
Haniiltonians are hermite (if we neglect decay widths of boundstates) and the boundstate spectrum is invariant 
under this transformation. For our effective Hamiltonians, however, the transformation is not unitary (the 
charge is not hermite), since the Hamiltonians are dependent on the cm. energy P° of the system. 



to determine all the coefficients. The problems are also related to the use of the equation of 
motion to simplify the Lagrangian since the on-shell condition is the equation of motion for 
an asymptotic field. For comprehensiveness we prove in an appendix that it is justified to use 
the equation of motion to simplify the local NRQED/NRQCD Lagrangian and also to simplify 
local current operators; to the best of our knowledge such a proof has never been provided 
explicitly, although similar proofs for other effective field theories have been given |2^, |29| and 
the claim itself is widely accepted already. 

Below we will use the language of QCD consistently; nevertheless all of our arguments 
hold also for the QED boundstates. Throughout the paper we neglect non-perturbative effects 
(those effects which are typically parametrized by Aqcd) and restrict ourselves to arguments 
which can be understood from a summation of perturbation series in as to all orders. 

The organization of this paper is as follow. After reviewing general aspects of gauge de- 
pendence of the conventional relativistic qq boundstate formalism (Sec. H), we summarize 
characteristic properties of the nonrelativistic boundstates from the viewpoint of the leading 
Coulomb singularities: their gauge independence and some nontrivial features are explained 
(Sec. Q). Then we define the effective Hamiltonian H^e for a qq system and investigate its 
gauge dependence as well as gauge dependences of the spectrum and wave functions of the 
boundstates using the BRST symmetry, within the framework of perturbative expansions in 
1/c; in particular we define the transformation charge Q of the effective Hamiltonian (Sec. ^, 
App. y). We clarify possible problems in the determination of i^eff if one uses the on-shell 
matching procedure (Sec. ^. For illustration, we present a calculation of the charge Q at 
(!?(l/c), corresponding to an infinitesimal gauge transformation from the Coulomb gauge; we 
also examine gauge cancellations in the decay amplitude of a qq boundstate into q'q"W^W~ 
at the same order (Sec. H). Conclusion and discussion are given in Sec. |^. In App. ^ we give 
a proof to justify the use of the equation of motion to simplify a local NRQCD Lagrangian. 
Some detailed discussions are given in App. |B| and |C[ 

2 Gauge Dependence of the Relativistic Boundstate For- 
malism: General Aspects 

We consider the BRST invariant QCD Lagrangian 

£ = -^ tr [G^^^G,,,] + Y, H^Pi^) - ^f] i^f + -^GF+FP, (1) 

/ 

where generally the sum of gauge-fixing and ghost terms can be written in a BRST exact form 
as 

Cgf+fp = {iQB,tT[cF]} (2) 

with the BRST charge Qb and an arbitrary gauge-fixing function F = F{A, ip, ip, c, c, B) |5^.[| 



* In our convention, the BRST transformation is defined as {iQs,')/'} = Sb^ — igcip, Sb^P — ^sV'c, (^s^^ 
-DpC, 5bC = igc^, Sbc = iB, and 5bB = 0, where B is the Nakanishi-Laudrup auxihary field. 



P/2 + P P/2 + Q 




P/2-P P/2-q 



Figure 1: The four-point function G{pqP). P/2 it q denotes the four- momentum of incoming 
quark/antiquark; P/2 it p denotes the four-momentum of outgoing quark/antiquark. 



Define a four-point function (Fig. P as 

G{xi,X2,X3,Xi) = {0\T^{xi)'tp{x2)^{x3)^{x4)\0) (3) 

_ f d^p dU d^p c(vaP) 



X exp 



-iy ■ (xi + X2 - X3 - X4) - ip ■ {xi - X2) + iq ■ (x^ - X4) 



•(4) 



Here and hereafter we restrict our discussions to a quark/antiquark of some specific flavor and 
omit tlie flavor index /. Tlie field operators and states are those of the Heisenberg picture. A 
quark- ant iquark boundstate contributes a pole to the Green function G{pqP). In the vicinity 

y]P) (with quantum number z/, mass M^, and 



of a pole corresponding to a boundstate 
momentum P), the Green function takes a form |]3T| 



where 



G{pqP) = ^-— po_,, ' .^ + (regular as P" ^ u^^p), (5) 

^"^!/,P -' ^V,P "•" ^^ 



(0 I T ^{x,)^{x2) I ^; ^ ) = e-^Kp^°-^^-^) y ^ e-^P- x.,p(p), (6) 

X.,p(p) = X.,p(p)^(7°®7°), (7) 



X = i(Xi+X2), X = Xi-X2, CU^p=yP2 + M2. (8) 

In this paper we assume that the decay width of a boundstate is infinitesimally small except 
where it is stated otherwise. 

An infinitesimal deformation of the gauge-fixing function, F ^> F + SF, induces a change 
of the Lagrangian 

f d^xSC = {iQB,SO}, SO= f d^xti[cSF]. (9) 

Accordingly the Green function changes as 

5G{xi,X2,X3,X4) = -{0\T{QB,5O}ij{xi)i^{x2Mxs)ij{x,)\0) (10) 

= -(0|T<5O[Qs,^(xiMx2Mx3)^(x4)]|0). (11) 



Suppose we are interested in the boundstates which can be created from the vacuum via 
gauge invariant operators, e.g. ip{^)4'{^)j ip{x)pip{x), ip{x)'^^ip{x), etc. For example, in the 
above Green function we may set xi = X2 = x, xs = X4 = y and contract color indices to make 
color singlet operators ^(a;)-?/'(x), il){y)il){y) . It then follows from [QbiiP'4'\ = that 

6G{x,x,y,y) =0. (12) 

color singlet 

Comparing this with eq. (§), one sees that both the mass and the residue of any boundstate 
which couples to the operator tptl) are invariant: 

5M^ = (13) 

and 

(5(0|V^(x)^(x) z/;P\ = 0, ie.5fAx.pip) =0. (14) 

/ J color singlet 

Since both | 0) and iIj{x)iIj{x) are BRST invariant, it implies that the boundstate satisfies the 
physical state condition 



Q 



B 



U-P)=0. (15) 



Note, however, that generally the boundstate wave function Xvpip) depends on the gauge- fixing 
condition: 



U-P\ 



^X.Ap) = -F.T. (0|T,5O[gB,V^(x)V^(i/)] 

= -F.T. (0|T(7 5O[c(x)V'(x)?(i/) + V^(x)^(l/)c(y)] z/; P ) , (16) 



where F.T. stands for an appropriate Fourier transform. 

Any physical amplitude (/;out|i;in) which involves the quark-antiquark boundstate con- 
tributions includes the above Green function G{pqP) as a part of it. Since the initial and final 
states satisfy the physical state conditions Qb U; in ) = Qb | /; in ) = and the theory is BRST 
invariant, the amplitude is gauge independent. Hence, the boundstate poles included in the 
amplitude are also gauge independent. An interesting question is whether the Green function 
G{pqP) includes any unphysical pole, which does not contribute to the physical amplitude, 
close to or degenerate with one of the physical boundstate poles. [] As for nonrelativistic quark- 
antiquark boundstates the answer is no, as will be shown in Section ^ and in Appendix y. 

3 Nonrelativistic Boundstates: Leading Coulomb Sin- 
gularities 

It is well-known that, in describing a system of a nonrelativistic color-singlet quark-antiquark 
(qq) pair, naive perturbation theory breaks down due to formation of boundstates ||32|, |33| . 

^ A typical example is the i?^-gauge for electroweak interaction where an unphysical pole (fc^ — ^M^ + ie)~^ 
is included in the gauge boson propagator. 
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Figure 2: The ladder diagrams for the process 7* -^ qq. The diagram where n uncrossed gluons are 
exchanged has a behavior ~ (as/P)" near threshold. 



Intuitively, this is because the slow q and q are trapped by the attractive force mediated by 
exchange of gluons and multiple exchange of gluons between the pair becomes significant. We 
review this property in a production process of a qq pair. 

Consider the amplitude of a virtual photon decaying into q and q, 7* -^ qq, just above the 
threshold. As we will see below, the ladder diagram for this process with n gluon exchanges 
has a behavior ~ (a^//?)"', see Fig. 0. Here, (3 is the velocity of g or g in the cm. frame. 



P 



1- 



Am'^ 



:i7) 



Hence, the contribution of the n-th ladder diagram will not be small even for a large n if 
/3 ^ as- That is, higher order terms in as remain unsuppressed in the threshold region. 
These singularities in (3 which appear in this specific kinematical configuration is known as 
the "Coulomb singularities" or "threshold singularities". The singularities arise because, for 
a particular assignment of the loop momenta, all the internal particles can simultaneously 
become almost on- shell as /? — ;> 0. 

The appearance of the factor (a^//?)" may be seen as follows. First, consider the one-loop 
diagram. Its imaginary part can be estimated using the Cutkosky rule (cut-diagram method), 
see Fig. ^ Namely, the imaginary part is given by the phase space integration of the product 
of tree diagrams. The intermediate qq phase space is proportional to /5 as 



d^2{qq) 






d cos 9, 



[li 



where 9 is the angle between the momenta of the intermediate and final quarks in the cm. 
frame. The qq scattering diagram with a t-channel gluon exchange contributes a factor ~ as/ P"^ 
since the gluon propagator is proportional to 1//?^; the propagator denominator is given by 



|^|2 = _^(l_COS^), 



(19) 




d&(^^) 



p 




1 X VP' 



Vp 



Figure 3: The Cutkosky rule for evaluating the imaginary part of the 1-loop diagram. The factors 
in as and [3 are shown explicitly. 




d $.(^S 




P X a^p X 



a/p2 _ (a/p)2 



^/^C. 



Figure 4: The cut-diagram method for evaluating the singularities of the higher order ladder 
diagrams. The factors in as and /5 are shown explicitly. 



where k denotes the gluon momentum. Thus, we see that the imaginary part of the one-loop 
diagram has the behavior ~ (5-asl i? = ots/ 13- Analyticity implies that the real part of the one- 
loop diagram has the same structure ~ as/ 13. By repeatedly using the cut-diagram method, 
one can show by induction that the imaginary part of the ladder diagram with n uncrossed 
gluons behaves as ~ (as/P)"", see Fig. §. 

Alternatively this fact can be shown by a power counting method f^. The relevant loop 
momenta in the loop integrals are in the nonrelativistic regime: 



P 



m, p° 



m 



k^ 






kr^O{(3). 



(20) 



Here, p, p and k represent the internal momenta of g, q and the gluon, respectively, in the cm. 
frame. For such configurations, g(g) and gluon propagators are counted as ~ l//9^, and the 
measure for each loop integration d'^k/{2'KY as ~ /^^.Q 

Thus, the ladder diagrams contain the leading singularities ~ {as/ (3)"'. Other diagrams (in 
particular crossed gluon diagrams) do not possess the leading singularities but only non- leading 
singularities ~ a^+7/5" (/>!)• 



* In counting the powers of /3 of a loop integral, the singularity of the integrand will increase if one assigns 
a large power of /? to the momentum in the propagators, but the integration measure is more suppressed. The 
optimal assignment of the order in (3 to each internal momentum must be sought to identify the most singular 
part of the integral. This procedure leads to Eq. (KG). 
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+ 




Figure 5: The self-consistent equation satisfied by the leading singularities of the qq'y vertex T^. 
One should take only the leading part ~ {as/ (i)^ on both sides of the equation. 



As the higher order terms in as cannot be neglected near threshold, we are led to sum up 
the leading Coulomb singularities. Let us first discuss gauge dependence of the amplitude when 
this summation is performed, in particular because only the ladder diagrams are included. The 
exact amplitude for the process e+e~ -^ qq near threshold can be expanded in terms of as and 

(3 as 



oo 

A^(-^"")(as,/5) = J2 Cn{as/I3T + (non-leading terms). 



(21) 



n=0 



The full amplitude A^(-^"") is gauge independent, so must be the each coefficient c^. Because 
only the ladder diagrams possess this type of singularities, the most singular part of the ladder 
diagrams has to be gauge independent. 

To see this explicitly, we examine gauge dependence of the gluon propagator. In a general 
covariant gauge, the qq -^ qq scattering amplitude in the threshold region is given by 



up^Ui 



—I 



kP' + ie 



Qtiu 



(1-0^ 



—I 



va Vf = UfYUi 



-k"^ + ie 



Vi'J^Vf X 



1 + OiP) 



(22) 



where the subscripts i and / stand for the initial and final state, respectively. We have used 
the fact that the space components of the currents, Hf'-f'^Ui and Vi'-f'^Vf, are order (3 in the 
cm. frame.[| Note that the leading part of the gluon propagator is identical with the Coulomb 
propagator in the Coulomb gauge. Eq. ( ^2]) also holds for the momenta ( pOf ) if we note that 
the off-shell q and q wave functions are given by 



: + m = m(l + 7°)+0(/3), 
^ + m = m(l-7°) + C(/5). 



(23) 
(24) 



Thus, gauge independence of c„'s is ensured by gauge independence of the leading part ~ l//3^ 
of the gluon propagator in eq. (^) . 

Let us denote by T^ the sum of the leading singularities of the vertex 7* -^ qq. It satisfies 
a self-consistent equation as depicted in Fig. ||. Retaining only the leading part ~ (as/P)"' on 

^ Dirac representation of the 7-niatrices is most useful in power countings, where 7° is diagonal and 7"s 
are off-diagonal. The quark spinor wave function has the upper two components of 0(1) and the lower two 
components suppressed by /?, vice versa for the antiquark. 



both sides of the equation, one obtains the vertex T^ as 



r" = - ( ^^r^-^ ) (E - p^m) G{p; E), (25) 



where E = \/s — 2m is the energy measured from the threshold. G{p; E) is the Green function 
of the nonrelativistic Schrodinger equation with the Coulomb potential: 

— - Cp^]-iE + te)] G(f; E) = S%r), (26) 

m r J J 

G{p; E)= f d^re~'P-''Gif; E), (27) 



where Cp = 4/3 is the color factor. The analytic expression of G{f; E) is given in terms of 
the hypergeometric function and includes the boundstate spectrum below threshold, E < Q. 
Alternatively, we may write 

Gi^-,E) = -t^^m^, (28) 

where 0n(p) and ipni'f^ are the Coulomb wave functions in momentum space and coordinate 
space, respectively. Here, n includes the boundstates with En = — (Ci7'as)^m/4n^ and the 
continuum states with £"„ > 0.[| 

At this stage, we see a nontrivial consequence of the summation to all orders in as- At 
any order of the perturbative expansion in as, the amplitude for 7* —>■ qq is zero below 
threshold, E < 0. For example, the absorptive part of a quark loop contribution to the 
vacuum polarization function 

(P V - P^P") Im Ug{P^) = I d^x e*^-(^-2') Im (0 | Ti){x)-^>'^{x)i){y)Yi'{y) I ) (29) 
vanishes below threshold. After summation of the leading singularities, however, it is given in 



terms of the Green function at the origin [32 



Imn,(.) = -^ImG(f = 0;E = v^-2m) = ^ ^ \i,M? ^E - E^), (30) 

n 

which in fact diverges at the positions of boundstates, E = En < 0. This discrepancy before 
and after the summation can be traced back to the fact that the limit e — i> in the propagator 
denominators does not commute with the summation to infinite orders in as- Namely, if we 
pursue the perturbative calculations with a finite e > 0, the absorptive part Imng(s) remains 
nonzero below threshold at each order. After the summation, constructive interference effects 
result in a drastic magnification of the amplitude ~ 1/e at E = En- 

In order to reach below the threshold for the process e"*"e~ —>■ qq, we need to include a 
subsequent decay process, e.g. q and q decaying into lighter quarks, or qq annihilating into 



^ To see that {E — p"^ /m)G{p\ E) is a function of as/ (i, one should identify E -^ mj3'^ and |p| -^ mfi at 
leading order. 
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multiple gluons, etc. Then the corresponding amplitude is nonzero and gauge independent 
both above and below the threshold. Summation of the leading singularities can be performed 
in the same way as above and leads to the same vertex V^, except that in this case the quark 
momentum \p\ needs not equal \/s/4: — m? (as required for an on-shell quark) as long as it is 
in the nonrelativistic region. | 

Below we discuss gauge dependences of the spectrum and the wave functions of the nonrel- 
ativistic qq boundstates within the framework of their calculations in perturbative expansions. 
An appropriate expansion parameter of this problem is 1/c, inverse of the speed of light, when 
c is restored as a dimensionful parameter [^. In this case both as = g'^/4:TThc and (3 = v/c are 
C(l/c) quantities.]^ Therefore the sum of the leading singularities (as/P)" is counted as C(l). 
Perturbative corrections to the boundstate wave function are given as a double expansion in 
as and /3, e.g. 0{l/c) corrections include a^+V/^" = asias/PT = I3{as/I3Y^^. Note that the 
parameter (3 is guaranteed to be small if as is small, since we are interested in the summation 
of the leading singularities only in the kinematical region where naive perturbation theory 
breaks down (/3 ^ as). The boundstate mass is given as a power series in as since they are 
independent of /3.| Throughout this paper we set c = 1 in our formulas in order to maintain 
simplicity of expressions; one may easily count the power of 1/c by counting the powers of as 
and (3. 

4 The Effective Hamiltonian for a qq System 

In this section we discuss gauge dependences of the spectrum and the wave functions of the 
nonrelativistic boundstates in a general framework. 

Let us introduce an effective Hamiltonian for a color-singlet qq system as follows. First we 
define a Green function for a qq pair in the cm. frame as 

g{p, q; A, A, A', A'; P°) = {p, -p, A, A I -— . I g, -g. A', A' ) , (31) 

where H denotes the full QCD Hamiltonian (including the gauge-fixing and ghost terms); 
I p, — p. A, A ) is an eigenstate of the free Hamiltonian Hq = H\ac;^o and represents a color- 
singlet two-body state composed of a free quark-antiquark pair: 

|p,-p,A,A> = a;^/ ,|0)f,,, ^o|0)f,,, = 0. (32) 

'^ color singlet 

Here, a^ {b'^) denotes the creation operator of a free quark (antiquark); p (— p) and A (A) 
denote the three momentum and the spin index of q (g) in the cm. frame, respectively. P'^ 

^ The nonzero decay width of the boundstate r„ renders the (5-function in eq. ( |30| ) to the Breit-Wigner 
distribution 

r„/2 



TrS{E - E„) 



[E-E^Y+Tl/A 



^ Here, (3 symbohzes both \p\/mc and yjE/rru? for a nonrelativistic off-shell quark/ antiquark. 
" In addition to powers of as and /3, there appear also powers of logag and log/3 in these perturbation 
series. 
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represents symbolically the cm. energy of the qq system, but we take the three energies, P", 
2iy^2^rm^ an(j 1\/q'^ + vn^i not necessarily equal to one another. Note that the above two- 
body state is not a physical state, Qb \ V-, ^V-, -^j ^ ) 7^ 0, which stems from the fact that ifo 
is not BRST invariant. Then we define an effective Hamiltonian which operates only on the 
subspace spanned by the two-body states such that it generates the same Green function: 

Qip, q- A, A, A', A'; P°) = (p, -p. A, A | ^Ipo^^- I ^^ -•?' ^'' ^' > • (^3) 

Namely, the effective Hamiltonian (a quantum mechanical operator) is defined by 

i7eff(n=P°-^-'(P°), (34) 

where Q'^ip, q; A, A, A', A'; P°) = (p, —p, A, A | Q^^{P^) | q, —q, A', A' ) is the inverse of the Green 
function restricted to the two-body subspace [take the inverse of Q{p, q; A, A, A', A'; P°) consid- 
ering it to be a matrix with indices [p, A, A) and {q, A', A')]. For analyzing the nonrelativistic 
qq boundstates, one first calculates the effective Hamiltonian in a series expansion in 1/c, then 
uses ordinary perturbation theory in quantum mechanics for calculating the spectrum and the 
wave functions of the boundstates in perturbative expansions in 1/c. As we have seen in the 
previous section, the leading order Hamiltonian is given byQ 

i/ir = 2m + ^-C^^. (35) 

m r 

Let us briefly explain the background why we introduced the Green function, eq. (^). Sup- 
pose we consider contributions from a qq boundstate to some physical process. In a calculation 
of the corresponding amplitude using time-ordered (or "old-fashioned") perturbation theory, 
the above Green function always appears as a part of that calculation. This is parallel to the 
fact that the four-point function eq. (D appears as a part of the calculation of the same ampli- 
tude using the (Lorentz covariant) Feynman rules. Time-ordered perturbation theory is often 
more suited for calculations of nonrelativistic processes because additional quark-antiquark 
pair productions are suppressed by powers of 1/c. 

The rules for time-ordered perturbation theory are ||5^: draw time-ordered diagrams (e.g. 



time flows from right to left), assign a matrix element {i\Va\j) at the time of each vertex, 
and assign a propagator 1/(P° — Ei + it) to an interval between two adjacent vertices. Here, 
Va is an interaction term, H = Hq + N, Kkj -P" is the total energy of the system; | i ) and Ei 

a 

denote the eigenstate and the eigenvalue of the free Hamiltonian, respectively, Ho\i) = Ei\i). 
Then we sum over all the intermediate states, where in general the energy is not conserved, 
E, + po. 

Although there are many ways to derive the rules (see Appendix |BD, simple correspondences 
to the ordinary Feynman rules may be seen by integrating over the time components of loop 
momenta and over the time components of external particles' momenta of a Feynman diagram 
for an unamputated Green function. In Coulomb gauge, decomposing the quark and transverse 



* Presently the QCD effective Hamiltonian is known up to 0(l/c^) in Coulomb gauge, see e.g. M, W\ 
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-q-k 



Figure 6: A time-ordered diagram which contributes to Q. The dashed hne represents the instanta- 
neous Coulomb gluon; the wavy hne represents the transverse gluon. 



gluon propagators as 



A+(p)7° + 



cjp-i (m — p • 7)7° 



A„(p)7° 



UJ. 



Vp + 






rn?, A-i-(p) 



2ujf, 



2\k\ 



rV'injj 



«j 



1 



1 



\k\^^^k^-\k\+te k^ + \k\-ie 



(36) 
(37) 
(38) 



and using the Cauchy theorem, every wave function becomes on-shell [e.g. "^^ u{p, \)u^{p, A) = 
A_|_(p)] whereas the energy conservation is violated. The ghost propagator can be handled 
similarly to the transverse gluon propagator. Integrating the Coulomb propagator is trivial 
because it is independent of the energy of the gluon. By way of example, the diagram in Fig. |^, 
which contributes to the Green function Q, is given by 



{Cp ■ 4710:5') 



d^k 



(P0-2o;j7+ze)(P0-2a;g-+ie) 7 {2nf v/^' ^ ' ' v/- ; vy, ; 

^ - 1 / PA-A —1 

1 
X . 

(PO - u^ - uj^^.^ - uj^ - uj^^^ + ie)(PO - uj^ - uj^^^ -\k\^ie) 



(39) 



We return to the discussion of the Green function Q and the effective Hamiltonian ifefr- If 
we vary the gauge-fixing function, the QCD Hamiltonian changes as^ 



bO= \ d^xix[cbF\ 



H-.H-{tQB,SO}, 
and the corresponding change of the Green function is given by 

1 r.^ .^. 1 



SQip, q- A, A, A', A'; P') = - (p, -p. A, A 



Wb,SO} 



pO-H + te 



(40) 



g,-g,A',A'> 



pO-H + ie 

^ Here, we assume that SF is independent of SqA^, Oqc, etc.; otherwise the change of the Hamiltonian takes 
a different form. See Section o for a more general case. 
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(41) 



- (P, -P, A, A I IQb -^q 7^-— SO — . I q, -q, A', A' ) 



- (P^ -P^ ^' ^ I po_ff + ,, ^O ^o_^^,, ^Qb I g, -g, A', A' > . (42) 

Since Qb \ P, —p, A, A ) ^ 0, generally 5Q ^ 0, so the corresponding effective Hamiltonian 
also changes, SHes 7^ 0. As we use the effective Hamiltonian to calculate systematically the 
boundstate spectrum and the wave functions in perturbative expansions in 1/c, we would like 
to see how they depend on our choice of gauge. 

The mass M,^ of a boundstate is given as the position of a pole of the Green function Q. 
Equivalently, it is calculated from ifefr(-P°) by solving 

[M,-H,s{M,)]\iy;eS)=0. (43) 

In the discussion to follow, we consider only those boundstates which appear already at leading 
order for | z/; eff ).0 From the definition of M,y above, one may evaluate its deviation when the 
effective Hamiltonian is varied infinitesimally: 

c^ ^ ^;effUgefr(M,) \u;eS) 

" {u;eS\l-H'^^{M,)\u;eSy ^ ' 

In the numerator, the variation of the Hamiltonian can be written as 

= [P° - ifefr(P°)] 5Q{P') [P' - H^siP')] (45) 

according to eq. (0). Eqs. ( p4|) and (^51) imply that 6Q should contain a double pole (P° — 
Mj, + ze)~^ in order to generate a nonzero shift of the mass 5M^ 7^ [^. 6Q contains, however, 
only a single pole (P° — M^, + ie)~^, since the state 

^O m !t ^ ■ 'Qb I P, -P, A, A > (46) 

in eq. (^^ does not include the boundstate pole. This follows from the physical state condition 
eq. (p!5|) . Also one can see explicitly by power countings of diagrams that at any order of 
1/c expansion the above state does not contain this boundstate pole; a proof is given in 
Appendix |C[ Thus, the boundstate mass is gauge independent in spite of the fact that the 
effective Hamiltonian is gauge dependent. 

In addition, the proof in Appendix y also shows that there is no unphysical state which 
contributes a pole to the Green function Q that is degenerate with or close to one of the poles 
of the physical boundstates of our interest. Stating more explicitly, there is no unphysical 
boundstate with a binding energy ~ a\m. 

■i- Note that at leading order all boundstates in the spectrum are the Coulomb boundstates which are 
physical states. (In particular, all these states can be created from the vacuum via gauge invariant operators.) 
It already suggests that to all orders of 1/c there are no unphysical boundstates in the spectrum of iJeff which 
are degenerate with these physical boundstates. 

14 



Now let us define a quantum mechanical operator Q{P^) (which operates only on the 
subspace of two-body states) by 

(p,-p,A,A|Q(P°)|g,-g,A',A'> 

^ ^ A", A" 

x^-i(g',g;A",A",A',A';P°). (47) 

Then Q{P^) does not include the boundstate poles (P° — My + ie)~^. Q can be interpreted as 
the generator of the transformation of gauge-fixing condition as seen from the relations 

5g = -iQg + iGQ\ (48) 

and 

5H,s = [i^eff(P°) - P°] ^Q{P') - tQ\P') [H,^{P') - P'] . (49) 

The last equation concisely represents the transformation of the effective Hamiltonian in a form 
which clearly shows the spectral invariance; cf. eq. (^If). One may easily see that the charge Q 
has following properties: in general Q is not Hermite, thus the transformation is non-unitary; 
Q vanishes at leading order of the 1/c expansion; beyond leading order, even at some specific 
order of 1/c the charge Q contains all orders of as due to the form of eq. (|47|). We will confirm 
these properties by explicit calculations in Section ^ 

Another method to verify gauge independence of the boundstate spectrum is as follows. The 
on-shell qq scattering amplitude can be calculated using the reduction formula of time-ordered 
perturbation theory, 

M,^^,q = lim (po - 2wp.)(p0 - 2uj^) g{p, q- A, A, A', A'; P°). (50) 

(See Appendix |B[) If this amplitude is analytically continued to an unphysical region, it 
exhibits a pole at the position of the boundstate, P^ = 2ujjj = 2uj^ — > My. If we expand the 
amplitude as a Laurent series at the pole 

p 

^<i^-^^^ = po_M ^^^ + (regular as P° ^ My), (51) 

and calculate the mass My in a perturbative series in 1/c, My should be gauge independent 
at each order of 1/c, since Jvtqq^qq is gauge independent at any order of perturbation series in 

as- 
Next we turn to the boundstate wave function, which is defined from a Laurent expansion 
of the Green function at P^ = My a^ 

g{p, q; A, A, A', A'; P°) = MpA,X)^MyA') ^ ^^^^^j^^ ^^ ^o _ ^^)^ (53) 
P^ — My + te 



§ At leading order of 1/c expansion, i' = [n, s, s) and 

eo)(p, A,A) = Mp)UX)UX), Mr^ - 2m- ^^^^f^, 

where ^s(A) — (A|s) is a two-component spin wave function. Expressions of ipy and My up to ©(l/c^) for the 
boundstates can be found in [|l6[ |l9| . 
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or equivalently, 

V9^(p, A, A) = {p, -p, A, A|z/; eff) (53) 

with a normalization condition 

{u;eS\l-H'^siM,)\iy;eS) = l. (54) 

Alternatively, from the original definition of Q, eq. (|31|), one may express 

V9^(p, A, A) = (p, -p, A, A z/;P = 0), (55) 



where 



jy\ 



P = 0) is the eigenstate of the full QCD Hamiltonian H; see Section ^. Then from 



eqs. (^) and ( ^8]) the variation of the wave function is given by 

6Mp,X,X) = -{p,-p,X,X\tQBjj W— -'^O u;P = o) (56) 

' iW,, — H + te I 



and 



b^^{p,\\) = -z[Q(M,)-(^,](p,A,A) 



A', A' 



when the gauge-fixing condition is varied. The last equation shows once again that Q can be 
interpreted as the transformation charge. 

Looking at eq. ( |56D one might think that it is possible to mix different gauges in calculations 
of decay amplitudes of the boundstate. Namely, one might take the wave function ipu{p, A, A) 
calculated in one gauge (e.g. Coulomb gauge) as the initial state wave function and calculate 
the rest of the decay amplitude in another gauge (e.g. Feynman gauge). Generally final states 
satisfy the physical state condition Qb \ f) = 0, so the above equation may suggest that such a 
calculation gives the correct result (the result of a consistent calculation in one specific gauge). 
This expectation, however, is wrong since the two-body states | p, —p, A, A ) do not span the 
complete Fock space. This fact will be verified explicitly in the second example in Section |^. 

5 Problems with the On-shell Matching Procedure 

In the definition of the effective Hamiltonian in terms of the full QCD Hamiltonian [eqs. ( |3TD 
and (P^j we kept the energies of the initial and final qq states different from P^. (It corre- 
sponds to off-shell initial and final states in the language of a Lorentz covariant formulation.) 
Accordingly the form of Hcs depends on our choice of gauge. In some literatures, however, 
the on-shell scattering amplitude eq. ( |5I1| ) is used instead of the off-shell Green function Q in 
order to determine a similar effective Hamiltonian. This leaves, in general, more freedom to 
the form of the effective Hamiltonian than what is due to gauge dependences. The difference 
is irrelevant when the Hamiltonian is applied to describe an on-shell qq system, whereas the 
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quark and antiquark inside a boundstate are generally off-shell. In this section we examine 
how the spectrum and the wave functions of the boundstates are affected when we employ the 
on-shell matching procedure to determine the effective Hamiltonian. 

First we consider a variation of the boundstate mass as we vary the effective Hamiltonian 
under the constraint that it gives the same on-shell scattering amplitude. As we have seen in 
eq. (^), SQ should include a double pole (P° — Mi^ + ie)~^ in order to generate a nonzero mass 
shift. We may try a simplest example: 

6g(p, q- A, A, A', A'; P") = (po _^mJ+ uf (2»)''''"(P - <!) ^xx'Syy, (68) 

Evidently it does not modify the on-shell amplitude (|50D, while it does generate a mass shift 
Ml, -^ Mj, + AMiy. In the calculation of the boundstate mass in a perturbative expansion 
in 1/c, if we add the above SHes to the effective Hamiltonian retaining terms up to some 
chosen order in 1/c, the mass is shifted up to the corresponding order. In fact one may find a 
variety of examples which can affect the boundstate mass while keeping the on-shell amplitude 
unchanged. Nevertheless we consider that it will not create a serious problem in practice, since 
we do not see any good reason why 6Hcs which has explicit pole structure (s) should mix in 
the determination of H^s- 

Next we consider the boundstate wave functions. Generally the wave function yj^ changes 
when SQ includes a single pole {P^ — M^ + ie)~^. For example, if we take 

sg = xg + gx, (60) 

i.e. 5H,s = [P°-ifeff(P°)]X + X[P°-i7eff(P°)], (61) 

the on-shell amplitude is not affected, where X is non-diagonal in momentum space and does 
not include the free particle poles (P° — 2uj^)'^, (P° — 2c<jg>)^^ On the other hand, the wave 
function varies as 

6ipi, = X ■ ip^. (62) 

In this case the variation is serious, since different calculations of a decay amplitude of a 
boundstate do not lead to a unique result if one uses different ip^^s connected by the above 
transformation as the intial state wave functions. 

One may think that the ambiguity related to the on-shell matching procedure to determine 
i^eff can be eliminated by directly matching all the relevant on-shell amplitudes to the pertur- 
bative expansion of the same amplitudes in as- This works at lower orders of 1/c expansions 
(in Coulomb gauge), but from the order 1/c'^ there appear contributions from the "ultra-soft 



gluons" which include all orders of as [§^] such that one should really consider the off-shell 
matching procedure seriously. 

We conclude, therefore, that the determination of the effective Hamiltonian Heg from the 
off-shell Green function g is favorable, and that the on-shell matching procedure can in general 
lead to incorrect calculations of the boundstate masses and the physical amplitudes involving 
boundstates. 
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6 Examples 



In this section we apply our formalism to two examples, where we study an infinitesimal gauge 
transformation from the Coulomb gauge. First example is a calculation of the transformation 
charge Q] in the second example we study gauge dependences of diagrams for a decay amplitude 
of a boundstate. 

Let us consider a class of gauge-fixing functions which interpolates the Coulomb gauge and 
the Feynman gauge. The gauge-fixing function is chosen as 

F = -2^ Q5 + d,A^ + ^^nV- a\ , (63) 



from which one obtains 



-^GF+FP — — tr 






1 „^ -\ 



(64) 



after integrating out the auxiliary field B. Here, .^ > is the gauge parameter: .^ ^ and 
(^ ^ 00 correspond to the Coulomb gauge and the Feynman gauge, respectively. The gluon 
propagator iD^^{k) is given by 

where a = 1 + ^m^/|A;p. Our formal arguments in the previous sections do not apply directly 
to this gauge-fixing condition since 6F includes OqA. Nevertheless we may obtain necessary 
rules for time-ordered perturbation theory easily via relations similar to eqs. (|36|)-(|38D.f^ For 
an infinitesimal change of the parameter ^ — >■ ^ + ^^, 

^0 = Jd'S^^tT[cnV-A]. (68) 



* A more natural choice of gauge-fixing function that interpolates the Coulomb gauge and the Feynman 
gauge would be 

F = -2i ( ifi + ^f,A^' + jV-A 

In this case, canonical quantization can be performed straightforwardly following the standard procedure pO| 
and all of our formal arguments apply directly. On the other hand, practical calculations are tediously compli- 
cated in this gauge due to the existence of a double pole (/c^ + ie)^'^ in the gluon propagator. For simplicity 
of practical calculations, we present the examples according to the gauge-fixing condition eq. (p3|). Another 
class of gauge-fixing conditions which interpolates these two gauges was introduced, for QED, in [Q, which 
corresponds to a class of nonlocal gauge-fixing functions. 
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Figure 7: The tree diagrams wliidi contribute to the charge Q for an infinitesimal transformation 
from the Coulomb gauge. The dotted line represents the ghost. The wavy line represents the gluon 
propagator zDjo; it is reduced effectively to an instantaneous propagator since the pole (/c^ + ie)~^ is 
cancelled by k"^ included in 50. 



6.1 The Charge Q at 0{l/c) 

First we calculate the transformation charge Q which generates an infinitesimal gauge trans- 
formation from the Coulomb gauge (^ = 0) at 0{l/c). For perturbative calculations it is 
convenient to rewrite eq. (^) as 

where H = Hq+V. P denotes the projection operator to the subspace spanned by the two-body 
states I p, —p, A, A ), and P = 1 — P. Time-ordered diagrams are obtained by expanding the 
above equation in terms of V and inserting the completeness relation in terms of the eigenstates 
of Hq. We may discard diagrams without cross talks between qq and ghost sectors, i.e. those 
diagrams which contain vacuum bubbles.[] The on-shell renormalization scheme is assumed for 
any value of ^, so we may neglect quark self-energy diagrams at 0{l/c). The BRST charge 
reads 

Qb= d^xg'ip\x)c{x)i!{x)^ , (70) 



where only the term which contributes up to 0{l/c) is shown. 

Simplest diagrams generated by eq. (|69|) are the tree diagrams shown in Fig. ^ The two 
diagrams give equal contributions, the sum of which is given by 

^ ' ' ' 0{as) \p - q]-^ P° - \p - q\ - ujj; - io^ + te 

Examining variations of the boundstate wave functions generated by this charge [eq. (0)], we 
see that two regions of the gluon-ghost momentum, 

soft : |p-g1~0(/5), 

ultra-soft : |p - g] ~ 0(/32), ^ ^ 

are relevant at (9(l/c).0 Existence of the ultra-soft region indicates that the diagrams with 

^ This corresponds to renormalizing the perturbative vacuum | 0)fj.j,j, appropriately in each gauge. 
'^ In power counting we consider d£,ni^ /\p— q\'^ ~ C'(l). 

19 



P-k 



Qb 




®50 

It-- 



®50 



-p 



-q-k 



Figure 8: The diagrams wiiicli contribute to the charge Q at 0{l/c) for an infinitesimal transforma- 
tion from the Coulomb gauge. The dashed line represents the Coulomb gluon of the Coulomb gauge. 
Other notations are same as in Fig. |^ 



multiple Couloinb-gluon exchange in ladder contribute also at 0{l/c). Indeed one may check 
that all the diagrams shown in Fig. |^ contribute to 5ipy at this order; the contributions come 
from the ultra-soft region of the gluon-ghost momentum, |fc| ~ 0{(3'^). This is also consistent 



with the result of Love 2t\. Hence, we find 



(p,-p,A,A|Q(P°)|g,-g,A',A'> 



0{l/c) 



i Cp 4770; c 



(2vr)3 1^1^ 



^('°Mp-|,g+|;A,A,A',A';PO-|A;|-^ , (73) 



where 



j(LO) 



(p, q; A, A, A', A'; P°) = (p, -p, A, A 



1 



PO-F 



(LO) 
cflf 



te 



g,-g,A',A'> 



(74) 



includes summation of the Coulomb ladders to all orders of as- The charge Q{P^) turns out 
to be anti-hermite at 0{l/c). We note that the above charge does not include any boundstate 
pole because of the integration over k. 

Alternatively it is possible to calculate the charge Q\o{i/c) by first evaluating 5Hes and 
then extracting Q via the relation eq. (^91). This procedure becomes cumbersome at higher 
orders of as because the number of gauge cancellations among diagrams increases. These 
gauge cancellations are automatically incorporated in the direct calculation of Q above by the 
BRST invariance of the full QCD Hamiltonian, [Qb-, H] = 0. 

6.2 A Decay Amplitude of a qq Boundstate at 0{l/c) 

Next we analyze infinitesimal gauge transformations of the diagrams for the decay process of 
the boundstate where q and q decay into lighter quarks via electroweak interaction. We analyze 
the infinitesimal transformation from the Coulomb gauge up to 0{l/c) as in the above example 
and see how the variation of the initial-state wave function eq. (0) gets cancelled in the total 
amphtude. The diagrams which contribute to this process up to 0(l/c) in Coulomb gauge 
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are shown in Fig. |^ |36[ .|^ When we vary the gauge- fixing function, additional diagrams which 
contribute to the 0{l/c) decay amphtude are shown in Fig. |TU|. Here, the double-wavy lines 
represent iSD^^, where 

etc. Diagrams (a) and (b) can be regarded as transformations of the initial-state wave function 
of the leading-order diagram (i). Conversely, the diagrams (c)-(e) cannot be regarded as such, 
since they do not contain a qq two-body state as an intermediate state. 

Using diagrammatic analysis, one may verify that the sum of all diagrams (a)-(e) van- 
ishes so that the total amplitude is indeed gauge independent. In fact, from diagrammatic 
manipulations as shown in Fig. |Tl] and also from similar manipulations corresponding to the 
diagram (b), one can show that the sum of diagrams in (a) and (b) can be regarded as the 
leading order diagram Fig. P(i) with the initial state boundstate wave function </?£, replaced by 
its infinitesimal transformation eq. (|57|) . Rearrangement of diagrams may be performed, for 
instance, using the relation 

1 / 1 1 1 \ 2(5fm2 1 



1 birri^ ( 1 

|^|2(po _ 1^1 _ ^^ _ ^^^^. + ,e) ~W V^° - 2^P + ^^ ' ^°-H-+fc + ^^ 



(76) 



for manipulating the propagator eq. ([75|). On the other hand, from Fig. |12| we see that the 
sum of the diagrams in (c)-(e) exactly cancels the sum of (a) and (b). For details of the 
diagrammatic analyses, see Refs. [^, |2^ . 

According to the formal arguments in Section ^ we know how the initial-state wave function 
transforms and therefore we know the sum of the other diagrams (c)-(e) in order to ensure 
gauge independence of the total amplitude. This example demonstrates that the diagrammatic 
analyses are rather cumbersome since infinitely many diagrams contribute even at the lowest 
nontrivial order of the 1/c expansion. 

7 Conclusion and Discussion 

In this paper we analyzed gauge dependence of an effective Hamiltonian formalism that de- 
scribes the nonrelativistic quark-antiquark boundstates and discussed problems of the on-shell 
matching procedure within this formalism. The significance of our present work may be put 
as follows. 

We used the BRST symmetry, which is known to be a powerful tool to study QCD Green 
functions, to analyze the NRQCD boundstates. The arguments were supplemented by power 



§ It is understood that the boundstate wave function ipv includes 0(l/c) corrections. For simphcity, we 
neglect 0{as) corrections to the qc(W and qq"W vertices, which constitute gauge independent subsets by 
themselves and do not mix with the gauge transformation of tp,^. 
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Figure 9: The diagrams which contribute to the amphtude for a nonrelativistic boundstate decaying 
into q'q"W^W~ up to 0{l/c). We suppressed diagrams for 0{as) corrections to the qq'W and qq"W 
vertices. 
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|5qv = -iQ-qv 
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(c) 



(d) 



(e) 



Figure 10: The diagrams which are generated by the infinitesimal variation of the gauge- fixing 
condition from the Coulomb gauge. The double-wavy line represents the variation of the gluon 
propagator i dD^^ . 




= z 






-4q-<r. 



Figure 11: A diagrammatic method to show that the sum of diagrams in Fig. |T^(a) can be regarded 
as the leading order diagram Fig. ^(i) with the initial state boundstate wave function ^pi, replaced by 
a half of its infinitesimal transformation —iQ ■ (fuf^- 
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etc. 



(c)+(d) + (e) = 



X 





X 



Figure 12: Examples of the diagrammatic method to show that the sum of diagrams in Fig. p!o|(c) 
exactly cancels the sum of (a) and (b). 
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countings of singularities of relevant diagrams to make them more explicit and detailed. Gauge 
dependence of the NRQCD boundstate formalism is more complex than that of usual (naive) 
perturbation theory since we have to deal with an infinite number of diagrams at each order 
of the 1/c expansion, (e.g. An infinite number of diagrams contribute to Hcs at 0{l/c) in 
gauges other than the Coulomb gauge |^.) We defined ifgfr naturally in the context of time- 
ordered perturbation theory. Then we obtained the transformation charge Q of Hes, from 
which we could easily see gauge independence of the spectrum and obtain transformation of 
the boundstate wave functions. For an infinitesimal transformation from the Coulomb gauge, 
we calculated Q directly up to (9(l/c). Also we saw that, without resort to the BRST symmetry, 
cumbersome gauge cancellations among diagrams are necessary to show gauge independence 
of a decay amplitude of the boundstate. At higher orders of 1/c, diagrammatic analyses such 
as what we presented in the second example or those in Refs. PSI, P^ become quite intricate 



so that the arguments based on the BRST symmetry would become more important. 

Furthermore, we showed possibilities for incorrect calculations of amplitudes involving 
boundstates if one uses only the on-shell qq scattering amplitude to determine Hes- These 
problems do not occur if we determine i^eff from the off-shell Green function Q, or, if we use 
a local NRQCD Lagrangian consistently and determine its coefficients via proper matching 
procedure, e.g. as in lattice calculations [^, |3^. The latter procedure has a disadvantage that 
one should calculate a number of amplitudes to determine all the coefficients. 

Presently we still do not have at our disposal a completely systematic way to identify all 
the necessary contributions in computations of physical quantities of the NRQED/NRQCD 
boundstates at a given order of 1/c expansion. We believe that the formalism developed 
in this paper will provide useful cross checks in these computations. Now we know how a 
boundstate wave function or the Green function Q contained in an amplitude transforms. 
The transformation charge Q is process independent and depends only on the gauge-fixing 
condition, and it can be calculated directly in a perturbative expansion in 1/c. 

A possible application is to use the formalism to study gauge dependences of the diagrams 
involved in the calculation of the top quark momentum distribution in the ti threshold region 
at (9(l/c^). It is known that at leading order the top momentum distribution is proportional to 
the absolute square of the wave functions of (would-be) toponium boundstates in momentum 
space [^. As we saw in this paper, wave functions of boundstates are gauge dependent beyond 
leading order. In the second example of Section H, we verified that this gauge dependence is 
cancelled by that of the final-state interaction diagrams (ii)-(v) at 0{l/c). In other words, a 
boundstate wave function mixes with the final-state interaction diagrams by gauge transforma- 
tion. This shows that the present calculations of the top momentum distribution p^ are gauge 
dependent, i.e. they vary if we transform the gauge infinitesimally from the Coulomb gauge, 
since they do not include the final-state interaction diagrams. Also the example suggests how 
gauge cancellations should take place in the complete amplitude at (9(l/c^) which has not been 
obtained yet. 
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A Use of the Equation of Motion in a Local NRQCD 
Lagrangian 

In writing down a local NRQCD Lagrangian in terms of the nonrelativistic quark (ipg) , nonrel- 
ativistic antiquark (ipq), gluon (A^), ghost (c) and antighost (c) fields, in principle one writes 
down all possible local interactions consistent with the rotational and BRST symmetries. In 
addition, one may simplify the Lagrangian using the equation of motion, and it is often con- 
venient to eliminate all terms including Dq {n >2), where D^ = d^ — igA^{x) is the covariant 
derivative. After such a simplification, the Lagrangian takes a standard form: 

C = ijlix) UDo + C2^ + C4-^ + cp^B ■ a + c^-^ (d ■ E - E ■ D) 
+ cs^ia ■ (dxE-ExD) +■■ ■]^q{x) + {^q -^ ^q) 

- itr [G^'^G.A . (77) 

We suppressed the gauge-fixing and ghost terms. One should determine the (Wilson) coeffi- 
cients of local operators C2, C4, cp, etc. by matching various on-shell amplitudes to those of 
full QCD. Furthermore, in practical applications of the NRQCD formalism, we often evaluate 
the correlators involving the current operators composed of the nonrelativistic quark and/or 
antiquark fields. The equation of motion is also used to eliminate Dq from the current opera- 
tors, and the coefficients of local operators constituting the current operators are determined 
by matching the on-shell amplitudes to those of full QCD. 

In this appendix we prove that we may use the equation of motion appropriately in order 
to simplify the form of the Lagrangian. We also prove that in the evaluation of on-shell 
amplitudes involving current operators, the change of the Lagrangian can be compensated 
by local redefinitions of the current operators and that one can use the equation of motion 
to rewrite the current operators. It is understood that we regularize ultraviolet and infrared 
divergences using the dimensional regularization. 

Let us start from a general local Lagrangian C{ipq,ipq,A^,c,c) and add a local operator 
which vanishes by the equation of motion: 

C^C + ^l{N,M}ijq. (78) 

Here, the equation of motion for ipq is denoted by 

-^ = (Mijq) (x), S= [d'^xC, (79) 
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Figure 13: The diagrammatic representation of eq. (pS]). ® shows the position of the local operator 
A^. The pole position is not changed, while the Z factor changes. 



and N denotes a local operator with A^ = A^"!", e.g. N ifjq = iDoipg, D'^ipg, B ■ aipq, etc. A^ may 
include the gluon field but not the quark or antiquark field. For simplicity we do not change 
the antiquark sector of the Lagrangian in our argument. According to eq. (|78D , the two-point 
and four-point functions change as 



S{0\Tijg{x)ijUy)\0) = {0\Tijg{x)ijUy) 



I d''z^Pl{z){{N,M}^P,){z) 



6 {0\T^g{x)ijl{y)^,{x')^l{y')\0) 



'\j,U.,'\ 



{0\T^Pgix)^^Jliy)^gix')^Uy^ 



'\j,Uo,'\ 



t d^z^Pl{z)i{N,M}^Pg){, 



0), (80) 



0).(81) 



In order to rewrite the right-hand-side of eq. (|80D one may use the Schwinger-Dyson equation[ 

i [ d^zN{z,z) 



(0 I t{^Mx) [{Ni,g){y)]^ + i {Nijg){x) ijliy) + ^,(a;) i^Hy) 



+ ipq{x) ipUy) 



t d^zijl{z){{N,M}^q){z) 



}|0) = 0. 



(82) 



The third term of this equation vanishes within the dimensional regularization, since N{z, z) 
contains ^'^(O) and/or [d^S^{z)]z^o which give scaleless integrals (tadpoles). Hence, we have 

5(0|T^,(x)^t(^)|0) = _(o|T{^^,(a;)[(A^^,)(y)]t + i(A^^J(x)^t(^)||0). (83) 

This equation shows that the change of the Lagrangian does not affect the pole mass of the 
quark propagator, whereas the Z-factor (wave function renormalization constant) varies; see 
Fig. 0. 

Following similar steps, one can show that the variation of the four-point function is given 
by 

5 {0\T^q{x)^Pl{y)^P,{x')4{y')\0) 

= - (0 I T {i ^,(x) [{N^gM]^ + i (iV^,)(x) ^liy)}^P,ix') ^J-(y') | ) . (84) 



* In the path-integral formulation, this follows readily from 



vi^l 



d4{z) "' 



^,(x)^J(y)[(7V^,)(z)]te^^ = 



and a similar term with ip -i^ ip^ . 
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Figure 14: The diagrammatic representation of eq. (p4|). Tlie first two diagrams give rise to a 



wave function renormalization common to Fig. 13. The third diagram is one-particle irreducible with 
respect to the leg with (gi, hence it does not contribute to the on-shell amplitude. 



Thus, if we redefine the Z-factor according to eq. (^31), the on-shell amplitude of the quark- 
antiquark scattering remains the same; see Fig. |1^. Similarly the amplitudes where multiple 
gluons are attached to the quark-antiquark scattering can be shown to be invariant under the 
variation of the Lagrangian eq. ( [THD . 

Eqs. (^) and ( [5^ ) also show that, when evaluating correlators involving current operators, 
the change of the Lagrangian can be compensated by local redefinitions of the current operators. 
By way of example, for a current operator which creates and annihilates a quark-antiquark pair. 



T(x) 



V^J(x) 



cla' + cla' 



. D' 



12m2 



+ 



i/jUx) +h.c.. 



(85) 



the on-shell amplitude calculated from the correlator {0\T J^{x) tljq{y)tljq{z) |0) remains un- 
changed if we redefine the current as 



J\x) -* J\x) + 5J\x) 
5J\x) = [(iV^,)(x)]t 



cV* + dia' 



D^ 



I2m? 



+ 



iil{x) + h.c. 



(86) 
(87) 



Finally we show that we may use the equation of motion in order to rewrite the current 
operators. One may derive the Schwinger-Dyson equation^ 



(0 I T I [{rM^,){x)]'^4{x)^,{y) ^ljq{z) + ir\x, y) ^{x) ^,{z) | | ) = 



(88) 



where F* is a local operator and may include the gluon field but not the quark or antiquark field, 
e.g. F*(x,|/) = a''5^{x — y), D\x,y), etc. The second term does not contain the quark pole, 
hence it does not contribute to the on-shell amplitude. Thus, adding [(F*M-?/'q)(x)]'''^|(a;) +h.c. 
to the current operator J*(x) does not affect the on-shell amplitude. 



^ This follows from 



and integrating over z. 



V^} 



Hq{z) 



V^X, Z) 'tj)\{x) i^qiy) ll)q{z) C 



iS 
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B Time-ordered Perturbation Theory 

Here, we derive the rules for calculations of the on-shell quark-antiquark scattering amplitude 
in time-ordered (old-fashioned) perturbation theory. The S'-matrix element between the eigen- 
states of the free Hamiltonian defined in eq. (p2D with an infinite time separation (asymptotic 
states) is given by 



% = \im^{p,-p,X,X\e-'''^\q,-q,X',X') (89) 

= hm / ^ e-^^"^ {p, -p, A, A | \ \ g, -g, A', A' > . (90) 

In the integrand, we see the Green function ^(p, g; A, A, A', A'; P°) introduced in eq. (|3l|) . We 
expand the right-hand-side in V", where H = Hq + V, 

1 1 °° / 1 \" 

P'> - H + le ^ P'> - Ho + ie^ y P^ - Ho + ie) ' ^^^^ 

n=0 ' 

and insert the completeness relations in terms of the eigenstates of Ho- One readily sees that, 
at each order of the perturbative expansion, the free propagator poles (P° — 2u^+ ie)~^ and 
(P" — 200^-+ ie)~^ are attached at the both ends. Therefore, if we write 

^^' -^' ^' ^ I PO - g + ^e I ^' -^' ^ ' ^ ) = (po _ 2a;,+ .e)(PO - 2u,+ ^e) ^'^^ 



and set 2c<jp' = 2uj^ = y/s, we find 

'" - ."i'i/^^-""(p^3^)'-^-«(^")- (93) 

Suppose A4qq^qq{P'^) IS regular inside the integration contour. Then 



T^oo dP° 

5^/ A/. _ _ 



po^v^ 



^hm e-^v^^{>i;,^,,(v/^) - zT A^,,^,,(v/^)}. (95) 



The second term in the last line represents the dominant term as T — * oo. Taking into account 
the irregularity oi Aiqq^qq{P^), we obtain additional terms which are subleading in comparison 
to the second term of eq. (|95D as T — > oo. Thus, we obtain the reduction formula eq. (|50D 
as well as the rules for calculations of the scattering amplitude in time-ordered perturbation 
theory, as explained in Sec. ^.Q 

Following similar steps, one can show that in general the Green function Q appears as an 
intermediate matrix element when one evaluates a transition amplitude involving contributions 
from the quark-antiquark boundstates using time-ordered perturbation theory. 



* The phase factor e *vs^ always appears in a perturbative evaluation of Sfi. It is irrelevant if we are 
interested only in the absolute value \Sfi\. 
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Figure 15: Typical diagrams which have power countings /?"'^ (as/P)^ for n 3> uq. The ghost is 
connected with the uncrossed ladders of qq with a finite number [ ^ ©(no)] of lines. 



C Absence of Boundstate Poles in Eq. ( [46| ) 



We show that the state given by eq. p6|) cannot accomodate a pole which is degenerate with 
any of the quark-antiquark boundstate poles, (P° — M^ + ie)~^. We first note that Qb and 
SO have the ghost number +1 and —1, respectively. Suppose this state contains some of 
these boundstate poles. Then, the matrix element composed of this state should have a power 
counting in terms of as and /? as 

(p, -p, \,\\60 — — -—— IQb I g, -g, A', A' ) 

~ r Y, Cn {ois/^r X [l + 0(l/c)] (96) 

n 

for some no, since P^ — M^ = [1 — C|n(a5//9)^/4/^] mfS"^ at leading order. It is known that the 
diagrams which can have the leading power counting {as/ (3)^ are only the uncrossed ladder 
diagrams; see Section ^. Therefore the diagrams which can contribute to Z?""" {as/f3)^ for 
n ^ uq are only those diagrams where a ghost is connected to one of the uncrossed ladders 
of qq with a finite number [ ^ O{no)] of lines; see Fig. 15.] After integrating over the loop 



momenta, there remains no pole in the P*^-dependence of the sum of the diagrams, in the same 
way that a usual one- loop diagram does not exhibit a pole but rather contains branch point (s); 
cf. eq. (p. 

We may restate it differently. If a ghost and a nonrelativistic qq pair should constitute a 
boundstate, intuitively the sum of the ladder diagrams with multiple gluon exchanges between 
the ghost and qq pair may exhibit a boundstate pole. Since the coupling of ghost and gluon is 
suppressed by powers of P, the binding energy of the boundstate should scale differently from 
(have more powers of 0:5 than) the Coulomb binding energies (if the boundstate should exist 
at all). 



''' We discard the diagrams without cross talks between qq and ghost sectors; see Section | 
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